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OPTIMIZATION OF ADIABATIC MAGNETIC MIRROR FIELDS 
FOR CONTROLLED FUSION RESEARCH 
byJ. Reece Roth 
Lewis Research Center 

SUMMARY 

This report formulates guidelines for the design of magnetic mirror machines used 
in controlled fusion research apparatus. The following optimization problem is consid- 
ered: a plasma of density n Q and kinetic temperature V , which is isotropic in velocity 
space, is released at the midplane of a magnetic mirror machine. The dimensions, field 
strength, and mirror ratio of the magnetic mirrors are then adjusted in such a way as to 
adiabatically confine the largest possible number of particles per dollar of investment in 
the magnetic field coils, hi keeping with the current state of controlled fusion research, 
it is assumed that the numerical magnitudes of the plasma density and kinetic tempera- 
ture are dependent on the outcome of the experiment and are not independently adjustable 
parameters which the experimenter can vary at will. Results of a recent experimental 
investigation are used to provide criteria that will assure adiabatic confinement of the 
plasma ions. Scaling relations, derived for both conventional and superconducting coils, 
relate the cost of the coils to the coil dimensions and the magnetic field strength produced 
by the coil. Such a scaling law for superconducting solenoids was derived from data on 
98 superconducting coils that have been sold commercially. 

It is shown that under many conditions, including conventional water-cooled copper 
coils, there exists no optimum design, and the total number of adiabatically confined par- 
ticles is just proportional to the amount of capital one can invest in the magnetic field 
coils. For superconducting coils, however, there exist optimum values of the magnetic 
field strength, coil dimensions, and mirror ratio for which the number of adiabatically 
confined particles per dollar of investment in the coils will be a maximum. For super- 
conducting coils, the optimum mirror ratios lie between the high values 1. 00 < B m ^ n / 
B m „ v < 0. 95. This optimum mirror ratio is relatively insensitive to the scaling law pa- 
rameters for magnet cost, provided that superconducting coils are under discussion. 

This high optimum mirror ratio, and the other guidelines developed in this report, are 
not qualitatively dependent on the cost of superconducting coils at any given time. 
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INTRODUCTION 

Until the present time, it has not been customary to optimize systematically the mir- 
ror ratio, the maximum magnetic field strength, and the dimensions of a proposed mag- 
netic mirror machine intended for controlled fusion research. This state of affairs has 
come about for several reasons: 

(1) There has been no generally accepted figure of merit for mirror machines, the 
maximization or minimization of which is a desirable goal of the design process. 

(2) It has been widely realized that nonadiabatic particle losses, which result from 
an escape cone in velocity space larger than would otherwise exist (refs. 1 and 2), may 
be avoided by designing the apparatus with mirror ratios R = B . /B near unity, 
or with strong magnetic fields, or with small magnetic field gradients. Until recently, 
however, there has been no quantitative relation among these variables to tell the appara- 
tus designer exactly where nonadiabatic losses will start to occur. Figure 1 shows the 
adiabatic and nonadiabatic loss cones in velocity space. The particles whose velocity 
vectors lie between these two loss cones are those whose loss is attributable to nonadia- 
batic effects. 

(3) If an economic constraint is to be applied to the design problem, the cost of the 
magnetic field coils must be expressed as a function of the magnetic field strength and 



Figure 1. - Regions of adiabatic and nonadiabatic particle losses in velocity space. 
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^geometry. Such scaling laws have not been available for superconducting coils. 

In this report, a figure of merit is proposed, the maximization of which is consistent 
with the requirements of plasma stability, the minimization of particle losses, and the 
maximization of the product of plasma density and kinetic temperature. A constraint on 
this figure of merit is given by the results of recent experimental measurements cm non- 
adiabatic particle losses, which are summarized in a formula that predicts the mirror 
ratio, apparatus dimension, and field strength for which particles begin to be nonadiabati- 
cally lost (ref. 2). The second constraint is provided by a scaling law for the cost of 
superconducting and conventional magnets. 

The optimization process is treated as a calculus of variations problem in which the 
figure of merit is maximized, subject to the constraints imposed by the coil cost and the 
requirement of adiabatic particle confinement. This procedure yields general expres- 
sions for the optimum mirror ratio, magnetic field, and apparatus dimensions as func- 
tions of the fixed parameters of the constraint equations. The best available estimates of 
the values of these parameters are substituted into the general expressions, and the opti- 
mum values of mirror ratio, dimensions, and magnetic field are exhibited for mirror ma- 
chines designed under current conditions. The plasma density and temperature are left 
as floating variables to be determined by the characteristics or the outcome of the experi- 
ment. 

The author would like to acknowledge the cooperation of the following organizations 
and individuals, who provided unpublished data on the dimensions and cost of superconduct- 
ing magnetic field coils: The Westinghouse Electric Company; A. J. Donius, of the Mag- 
nion Corp. ; N. S. Freedman of the Radio Corporation of America; Z. J. J. Stekly of the 
AVCO Corporation; J. F. Howlett of the Cryonetics Corporation; and D. C. Freeman, Jr., 
of the Union Carbide Corporation. At the request cf certain of the organizations which 
supplied information, the names of the manufacturers are not associated with the specific 
coils listed in tables I to m. 

FIGURE OF MERIT FOR CONTROLLED FUSION RESEARCH APPARATUS 

The present state of controlled fusion research is such that a figure of merit can be 
only empirical at this time. However, a careful study of existing mirror machines and 
research objectives has led to the conclusion that the following factors should be taken 
into account in any proposed figure of merit: 

(1) The figure of merit should be proportional both to the plasma density n and to 
the plasma kinetic temperature V (kT = eV), since one of the basic objectives of con- 
trolled fusion research is to raise the energy density to a level high enough so that con- 
trolled fusion reactions take place. 
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(2) The ions that constitute the fuel of the thermonuclear reaction must, on the aver -4 
age, be confined long enough to react. The figure of merit should, therefore, be propor- 
tional to the average characteristic confinement time of an individual ion in cases where 
this time is equal to or less than the mean free time between fusion reactions at the exist- 
ing plasma density and temperature. 

(3) The figure of merit should be proportional to the volume of the plasma, since 
large dimensions promote adiabatic confinement (refs. 1 and 2) and reduce the severity 
of the flute instabilities that are driven by magnetic field gradients (ref. 3). In addition, 
large dimensions reduce the surface-to-volume ratio. This reduction makes plasma con- 
tamination less of a problem and makes the effect of inserting diagnostic probes relatively 
less disturbing to the bulk of the plasma. 

(4) Both theoretical (ref. 4) and experimental (ref. 5) investigations have shown that 
departures from a Maxwellian distribution along a radius vector in velocity space pro- 
mote the growth of plasma instabilities. Such departures should be avoided by choosing 
a suitable injection method or by promoting equilibration among the injected ions. The 
figure of merit may therefore be based on the assumption of a Maxwellian energy distri- 
bution. 

(5) Theoretical (ref. 6) and experimental (ref. 7) investigations have shown that an- 
gular anisotropies of the velocity distribution in velocity space promote the growth of 
microinstabilities. Such angular anisotropies are to be avoided for other reasons. The 
anisotropy characteristic of all magnetic bottle geometries is an escape cone through 
which particles are lost by collisional scattering, similar to that illustrated in figure 1. 
The smaller is the escape cone, the more nearly isotropic is the velocity distribution and 
the smaller will be the losses due to collisional scattering. The figure of merit should, 
therefore, contain some factor proportional to the angular isotropy of the velocity distri- 
bution in velocity space. 

The proposed figure of merit may be written as 

f r rty r 00 /~2ir o 

S2=/ / / / / 0 Vf(v)v sin 9 d 9 dcp dv di {/ dt (1) 

Jo Jo J o Jo Je Q 

where t is the average e -folding decay constant of the plasma ion density, p is the 
volume of the plasma, v, cp, and 9 are the spherical coordinates of velocity space, f(v) 
is the time -dependent velocity distribution function of the ions in velocity and configura- 
tion space, and V is the average particle energy in electron volts. (All symbols are 
defined in the appendix. ) 

In equation (l), it has been tacitly assumed that the only anisotropy in velocity space 
is an adiabatic, axisymmetric escape cone of half -angle 0 Q (fig. 1), whose axis is par- 
allel to the magnetic field lines in the plasma. 
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* It is assumed that the plasma is in the steady state and that the distribution function 
is Maxwellian along a radius vector in velocity space. The distribution function is then 


f(v, cp, 6, t) = f(v) = n / m \ ^ exp (- 6 - e 0 ( 2a ) 

\2ireV/ \ 2eV / 

and 

f(v, (p, 0 , t) = o e^e Q (2b) 

so that the loss cone is entirely free of particles. If equation (2) is substituted into equa- 
tion (1), the figure of merit may be written as 



exp 



sin 6 d0 dcp dv di^ 


(3) 


where the integration over confinement time has been performed. Integrating over the 
coordinates in velocity space gives 


n = 



nV cos 0 Q 


dip 


(4) 


It should be noted that the results of this analysis are not affected by the functional 
dependence of the figure of merit on plasma density and temperature. The energy density 
nV has been used in equation (4) for the sake of concreteness, thus making the figure of 
merit proportional to the total plasma energy confined in phase space during a mean con- 
finement time. The thermonuclear reaction rate n (a v) could also be used in equa- 
tion (4) without changing the conclusions of this study. Under these conditions, the figure 
of merit would then be equal to the total thermonuclear energy released by the plasma 
during a mean confinement time. 

E has been shown (ref. 8) that the loss cone angle in velocity space for a magnetic 

mirror with maximum field B is given by 

max 


cos 8 = 
o 



(5) 


where B is the local value of the magnetic field. K equation (5) is substituted into equa- 


5 


15^1 [SI EE3 


tion (4), the figure of merit may be writ - 4 


S' 


ten as 


\ 


si isi isi 




ft 


z o- 


"/"f 


B 


B 


di// 


( 6 ) 


max 
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The geometrical configuration of the 
plasma has been assumed to be that of the 
upper half of figure 2 in which a long, cy- 
lindrical, uniform -field center section is 
bounded on both ends by magnetic mirrors. 
As indicated in figure 2, the uniform field 
— * center section will be ignored for two rea- 
sons: the uniform field will assure adia- 
batic confinement of the particles in this 
region, and, in addition, the optimum 
placement of the center section coils can- 
not be found without a long and detailed design study for each geometry considered. The 
plasma volume assumed in the subsequent discussion is shown in the lower half of fig- 
ure 2. The limits of integration on equation (6) are, therefore, 


Figure 2. - Assumed geometry of magnetic field. 
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max 


where r^ is the mean plasma radius, and Z Q is the axial distance from the point where 

B is equal to B . to the point in the mirror throat where B is equal to B_„_. 

^ min max 

For axisymmetric mirror systems and for nonaxisymmetric magnetic wells of the 
Ioffe type, the variation of the total magnetic field along the axis is given approximately 
by 


B * -B 


max 


(1 + R m ) - (1 - R m )cos 


ttZ 


( 8 ) 


where R m is B .„ /B mov . 

m min' max 

After this definition of R m is substituted into equation (7) and 77 , a, and a Q are 
defined as 
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(9b) 


(9c) 


the figure of merit may be written as 


Vnrrjz /l - R X 1 / 2 9 / * 2,r 

q = ° 0 c ° ( 22. 1 k 2 


m JT 


Vl + cos tj dor d0 f drj 


( 10 ) 


where V Q is the ion temperature expressed in electron volts on the axis at the midplane 
of the magnetic bottle, and n Q is a fictitious density equal to the density of a distribution 
isotropic in velocity space located on the axis at the midplane of the magnetic bottle. 
Maximizing the figure of merit can then be accomplished by minimizing the escape cone 
angle. This would then increase the fraction of an initially isotropic velocity distribution 
which would be trapped between the mirrors. The actual density of particles is equal to 

n Q Vl - R m - The parameter K is the ratio of the mean plasma to the mean coil radius. 
A "geometry factor" is then defined by 



whose precise value depends on the detailed way in which particle density, particle kinetic 
temperature, and magnetic field depend on the dimensionless spatial coordinates a , 0', 
and tj. 

The geometry factor | is of order unity and, for the special case in which the energy 

density and the magnetic field are independent of position, £ has the value of 2^2 K 2 . 
This factor must be known so that the figure of merit of any two specific devices can be 
compared. It is sufficient to write the figure of merit as 
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ft = n V rr*,Z 

O O C 01 


a - R 


, 1/2 


m 


( 12 ) 


in the present study, since only its extremal value is of interest. The optimization of 
equation (12) requires a knowledge only of the functional dependence of the figure of merit 
on r c , Z Q , R m , and r and not the absolute value of ft. Equation (12) contains the de- 
sired product of the confinement time with the velocity space volume and physical space 

volume. The term Vl - R m in equation (12) is a measure of the size of the confinement 
region in velocity space. As discussed on page 4, this term should be made as large as 
possible to minimize losses due to collisional scattering into the escape cone and to re- 
duce the severity of the instabilities that are driven by angular anisotropies of the velocity 
distribution function in velocity space. 


DESIGN CONSTRAINTS ON CONTROLLED FUSION RESEARCH DEVICES 

Adiabatic Confinement 


The degree to which conditions are adiabatic may be measured by the value of the 
adiabatic parameter e, which is defined as 


= mv 

eB Z 
av o 


_L I 2mV 

J av Z o > e 


(13) 


where B ay is equal to (l/2)(B min + B max ). If e is small, the ratio of the radius of 

gyration to the apparatus dimensions is small, and conditions are adiabatic. If e is 

large, the conditions are nonadiabatic, and a magnetic mirror will not be fully effective 

in reflecting particles because of the enlarged escape cone. 

There are several reasons why e should be as large as possible, but not so large 

that nonadiabatic losses occur. For a given particle energy, an apparatus with large e 

will be less expensive than an apparatus with small e, since the cost is proportional to 

both size Z and magnetic field B . If, however, e is made too large, nonadiabatic 
u av 

losses will occur, and the apparatus will not confine particles effectively. The best eco- 
nomic compromise, therefore, is to make e as large as possible without making it so 
large that nonadiabatic losses occur. 

Stability considerations also imply a compromise value of e. It has been shown by 
Kuo, Murphy, Petravic, and Sweetman (ref. 9) that large values of e promote plasma 
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'stability by the mechanism of "finite Larmor radius stabilization”; therefore, e should 
be as large as possible. On the other hand, e should not be so large that nonadiabatic 
losses occur, since the additional anisotropy in velocity space associated with such losses 
will promote the anisotropy-driven instabilities discussed by Harris (ref. 6). The best 
compromise value of e, again, is to make e as large as possible without making it so 
large that nonadiabatic losses occur. 

A series of experiments has been performed recently to measure the critical value 
of e, at a given mirror ratio and radius in the magnetic field, above which nonadiabatic 
losses will occur (ref . 2). K was found that the critical value of e for a particular mir- 
ror ratio R m is given by 



(14) 


hi the axisymmetric geometry studied in the series of experiments in reference 2, the 
parameters A and e A for single interactions with a magnetic mirror were given by 

-0. 167a 

A = 0. 388 e ° (15a) 


and 


0.075a 

e A = 0. 348 e 0 


(15b) 


where a Q is the dimensionless radius defined by equation (9c). These experiments 

covered the range of parameters 0 ^ a Q ^ 1. 75, 0. 003 < e ^ 0. 3, and 

0.175 < R < 0.95. 
m 

In order to assure adiabatic confinement in axisymmetric magnetic mirror machines, 
therefore, the relation between the axial mirror ratio and the adiabatic parameter should 
be 


R - 
m 


-0. 167 a,. 
\0. 388e ° 


V0. 348 e 


0.075aj 


(16) 


where the inequality specifies the region of adiabatic confinement, and the equality repre- 
sents the experimentally derived boundary cf the nonadiabatic region. 

K was found from an auxiliary series of numerical computations (ref. 2) that, if all 
other conditions are held constant, the addition of multipolar (Ioffe) windings to an axi- 
symmetric mirror field will make the conditions within the apparatus volume less adia- 
batic. H; was also found that, if the apparatus is designed to be adiabatic at a given 
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radius, it will be adiabatic for all lesser radii. Equations (15) and (16) were obtained fo^ 
a single interaction of a particle with an axisymmetric magnetic mirror. In designing 
magnetic bottles that have superimposed multipolar windings and in which multiple reflec- 
tions occur, equation (16) must be treated as though it puts an upper bound on e and/or 
a lower bound on R , the exceeding of which will certainly result in nonadiabatic parti- 
cle losses at the radius under consideration. 


Scaling Laws for Cost of Superconducting Solenoids 

Most of the conventional steady -state magnet facilities currently in use in controlled 
fusion research consist of water-cooled copper solenoids, which absorb large amounts of 
direct current electrical power. Recent progress in the art of fabricating superconduct- 
ing solenoids has made it clear that such coils are not only substantially less expensive 
to build and operate than conventional magnet facilities, but they also offer many advan- 
tages of convenience and flexibility over conventional water -cooled copper coils (ref. 10). 
In reference 11, it is shown that superconducting solenoids are less expensive to build 
and operate than conventional water-cooled copper solenoids for all conditions of interest 
for controlled fusion applications (see figs. 3 and 4 of ref. 11). Progress since this ref- 
erence was published make superconducting solenoids even more attractive. 

In both = 25 kg (ref. 10) and B„„„ = 67 kg (private communication with Dr. C. 

max liioX 

Laver ick) superconducting magnet facilities, the capital costs of the superconducting mag- 
net facility were less than half that of an equivalent conventional magnet facility, and the 
steady -state running costs were no higher than those of a conventional facility. For these 
reasons, superconducting magnet facilities will be discussed more extensively than con- 
ventional magnets in this report. 

In order to maximize the figure of merit subject to an economic constraint, it is nec- 
essary to have a scaling law for the cost of a single coil as a function, for example, of 
B„„„, the maximum magnetic field on the coil axis, and d, the mean diameter of the 
coil. Information about the dimensions, the maximum magnetic field, and the cost of 
98 superconducting coils was obtained from commercial suppliers, and a scaling law was 
derived by obtaining a best -fitting curve to this data. 

The inner radius, outer radius, mean diameter, length of winding, coil volume, 
maximum magnetic field on the coil axis, and cost of 98 coils are listed in tables I, n, 
and III. The companies and individuals that supplied this data have been mentioned in the 
INTRODUCTION. Table I contains data from 6 coils wound with niobium -tin superconduc- 
tors, table II contains 85 single, plain coils wound with niobium zirconium, and table III 
lists 7 single niobium -zirconium coils enclosed in cryogenic Dewars. 

An attempt was made to fit the data in the tables to three candidate scaling laws, 
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( 17 ) 


$ = ?o + *l B max d 2 

* = *0 + *1 B ma* V ? < 18 > 

* - *0 ♦ *l B max + h* 2 < 19 > 

where B ra , v is the maximum magnetic field strength on the coil axis, and V is the 
volume of the superconducting windings. The data in tables I, n, and HI were used as in- 
put to a least -squares curve -fitting computer program that calculated the values of $q, 
and $2 for given values of y^ and The computer program also calculated the 

mean error 

N 

2 l $ k ' $ ak 

A. $ = M (20) 

1 N 


and the mean square error, 



( 21 ) 


TABLE I. - PROPERTIES OF HIGH -FIELD NIOBIUM-TIN COILS 


Coil 

Inner 

radius, 

r i> 

m 

Outer 

radius, 

r 2- 

m 

Mean 

diameter, 

d, 

m 

Coil 

volume, 

V 

m 3 

Axial 

length 

of 

winding, 

L, 

m 

Maximum 
field on 
axis, 

B max’ 

W/m 2 

Actual 

cost, 

1 

0. 0159 

0. 0858 

0. 102 

0.0022 

0. 098 

11.0 

30 000 

2 

.042 

.105 

.147 

.0035 

.121 

8.0 

100 000 

3 

.0762 

.242 

.318 

.0484 

.292 

15.0 

330 000 

4 

.254 

.406 

.66 

.0322 

.102 

4.0 

150 025 

5 

.254 

.406 

.66 

.0322 

.102 

4.0 

167 700 

6 

.254 

.406 

.66 

.0322 

.102 

4.0 

125 500 
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TABLE n. - PROPERTIES OF BARE NIOBIUM ZIRCONIUM COILS 


Coil 

Inner 

radius, 

r l> 

m 

Outer 

radius, 

r 2’ 

m 

Mean 

diameter, 

d, 

m 

Coil 

volume, 

V 

m 

Axial 

length 

of 

winding, 

L, 

m 

Maximum 
field on 
axis, 

B max’ 

W/m 2 

Actual 

cost, 

*A 

Coil 

Inner 

radius, 

r V 

m 

Outer 

radius, 

r 2' 

m 

Mean 
diameter , 
d, 
m 

Coil 

volume, 

V s. 

m 

Axial 

length 

of 

winding, 

L, 

m 

Maximum 
field on 
axis, 

B max’ 

W/m 2 

Actual 

cost, 

$ A 

1 

0.032 

0.089 

0.121 

0. 00145 

0.067 

3.2 

6 000 

46 

0.038 

0.060 

0.098 

0. 00148 

0.216 

2.50 

7 100 

2 

.076 

.146 

.222 

. 00294 

.060 

1.65 

7 980 

47 

.0064 

.035 

. 041 

. 00014 

.038 

3.00 

1 365 

3 

.065 

.158 

.223 

.00367 

.056 

1.5 

10 700 

48 

.0064 

.025 

.032 

.00011 

.060 



1 590 

4 

.138 

.214 

.352 

.00360 

.043 

.56 

18 500 

49 

.013 

.035 

. 048 

. 00023 

.070 



1 955 

5 

.026 

.051 

.077 

.00069 

. 114 

3.0 

4 000 

50 

.013 

.035 

.048 

. 00030 

.089 



2 745 

6 

.137 

.304 

.441 

.0706 

.305 

3.2 

125 000 

51 

.014 

.033 

.048 

.00044 

.155 

3.00 

2 850 

7 

.0064 

.029 

.035 

.00011 

.045 

3.0 

1 350 

52 

.025 

.051 

.076 

.00033 

.054 



3 980 

8 



.032 

.038 

.00023 

.076 

3.0 

1 950 

53 

.025 

.045 

.071 

. 00116 

.264 



7 095 

9 



.045 

.051 

.00047 

.076 

5.0 

3 950 

54 

.052 

.079 

.131 

. 00274 

.242 



18 450 

10 



.045 

.051 

.00074 

.121 

5.0 

4 500 

55 

.019 

.045 

.064 

. 00100 

.190 

3. 

50 

5 680 

11 

.0064 

.057 

.063 

. 00077 

.076 

6.0 

4 900 

56 

.013 

.048 

.061 

. 00107 

.159 

4. 

00 

6 150 

12 

.011 

.070 

.081 

.00238 

.159 

6.0 

12 500 

57 

.013 

.045 

.057 

.00082 

.143 

4. 

00 

4 965 

13 

.013 

.038 

.051 

.00041 

.102 

3.0 

2 700 

58 

.0064 

.043 

.049 

. 00027 

.048 

5. 

00 

2 560 

14 

.013 

.032 

.044 

.00034 

.127 

3.0 

3 300 

59 

.0064 

.038 

.044 

. 00034 

.076 

5. 

00 

2 830 

15 

.013 

.046 

.059 

.00078 

.127 

4.0 

5 650 

60 

.011 

.045 

.056 

. 00088 

.146 

5. 

00 

4 810 

16 

.013 

.076 

.089 

. 00191 

.108 

5.5 

6 700 

61 

.013 

.057 

.070 

.00080 

.083 

5. 

00 

4 535 

17 



.057 

.070 

.00099 

.102 

6.0 

7 100 

62 



.061 

.074 

. 00100 

.089 



5 085 

18 



.054 

.067 

.00132 

.153 

5.5 

9 500 

63 



.049 

.062 

.00077 

.108 



4 660 

19 



.054 

.067 

.00176 

.203 

5.5 

12 500 

64 



.051 

.064 

.00106 

.140 



5 665 

20 



.078 

.090 

. 00223 

.121 

8.0 

16 300 

65 



.051 

.064 

. 00135 

.178 



7 200 

21 

.013 

.076 

.089 

.00248 

. 140 

8.0 

18 890 

66 

.013 

.046 

.059 

.00097 

.159 

5. 

00 

6 770 

22 

.016 

.029 

.045 

.00017 

.095 

1.5 

1 100 

67 

.014 

.059 

.073 

.00064 

.064 



4 710 

23 



.067 

.083 

.00167 

.127 

6.0 

11 300 

68 



.059 

.073 

.00077 

.076 



4 895 

24 



.067 

.083 

.00234 

.178 

6.0 

14 500 

69 



.051 

.066 

.00063 

.083 



6 165 

25 



.067 

.083 

.00276 

.210 

6.0 

18 100 

70 



.054 

.068 

.00087 

.102 



6 165 

26 

.027 

.041 

.068 

.00016 

.051 

.20 

765 

71 

.014 

.051 

.065 

. 00105 

.140 

5. 

00 

6 950 

27 

.0254 

.039 

.065 

. 00021 

.076 

.50 

955 

72 

.018 

.052 

.069 

. 00113 

.152 



7 550 

28 

.095 

.108 

.203 

.00104 

.127 

.50 

3 010 

73 

.019 

.056 

. 075 

. 00103 

.121 



8 265 

29 

.0064 

.014 

.021 

.00002 

.045 

.80 

805 

74 

.025 

.067 

.092 

. 00196 

.165 



10 940 

30 

.014 

.029 

.043 

. 00013 

.070 

1.00 

955 

75 

.025 

.068 

.094 

. 00140 

.111 



8 250 

31 

.0222 

.038 

.060 

.00021 

.070 

1.00 

1 175 

76 

.013 

.056 

.068 

. 00134 

.146 

5. 

50 

7 715 

32 

.025 

.044 

.070 

.00106 

.254 

1.00 

2 390 

77 

.013 

.052 

.065 

. 00184 

.228 



9 485 

33 

.018 

.035 

.053 

.00043 

. 152 

1.50 

2 445 

78 

.014 

.060 

.075 

. 00167 

.155 



9 920 

34 

.018 

.035 

.053 

.00055 

. 194 

1.50 

2 745 

79 

.014 

.054 

.068 

.00087 

.102 



6 715 

35 

.020 

.033 

.053 

.00011 

.051 

1.50 

1 455 

80 

.018 

.060 

.078 

. 00253 

.241 

1 


12 090 

36 

.025 

.038 

.064 

.00026 

.102 

1.50 

2 560 

81 

.018 

.060 

.078 

. 00160 

.152 

5. 

50 

10 030 

37 

.025 

.038 

.064 

.00051 

.203 

1 

3 055 

82 

.0064 

.048 

.054 

.00058 

.084 

6. 

0 

5 290 

38 

.032 

.048 

.079 

.00079 

.200 


3 895 

83 

.019 

.067 

.086 

. 00204 

.159 


1 

11 530 

39 

.032 

.04 4 

.076 

.00069 

.228 


3 895 

84 

.025 

.073 

.098 

. 00252 

.171 



14 535 

40 

.013 

.032 

.044 

.00012 

.045 

2.00 

1 450 

85 

.025 

.083 

.108 

.00215 

.111 


k 

12 445 

41 

.013 

.032 

.044 

.00074 

.280 

2.00 

3 480 











42 

.025 

.041 

.066 

.00063 

.203 

2.00 

3 576 











43 

.038 

.060 

.098 

.00061 

.089 

2.30 

4 140 











44 

.027 

.048 

.075 

.00133 

.267 

2.50 

6 910 











45 

.029 

.046 

.075 

.00099 

.242 

2.50 

6 625 
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TABLE m. - PROPERTIES OF NIOBIUM -ZIRCONIUM COILS 


IN DEWARS 


Coil 

Inner 

radius, 

r u 

m 

Outer 

radius, 

r 2> 

m 

Mean 

diameter, 

d, 

m 

Coil 

volume, 

V s> 

m 3 

Axial 

length 

of 

winding, 

L, 

m 

Maximum 
field cm 
axis, 

B max’ 

W/m 2 

Actual 

cost, 

*A 

1 

0. 035 

0. 059 

0. 094 

0. 00105 

0.152 

2.0 

8 300 

2 

.038 

.076 

.115 

.00208 

.152 

5.0 

21 000 

3 

.029 

.070 

.099 

.00391 

.305 

5.0 

37 000 

4 

.070 

.134 

.204 

.00287 

.070 

2.5 

26 880 

5 

.083 

.147 

.229 

.00322 

.070 

2.5 

30 650 

6 

.095 

.159 

.254 

.00357 

.070 

2.5 

34 335 

7 

.095 

.159 

.254 

.00357 

.070 

2.5 

32 920 


TABLE IV. - BEST-FITTING PARAMETERS CALCULATED FROM DATA OF 

TABLES I TO m 


(a) For relation $ = $ Q + 


Table 

Number 

of 

coils 

Scaling law 
parameter 

Scaling law 
exponent 

Hoot 

mean 

square 

error, 

V 

Mean 

error, 

A x $ 

Mean 

coil 

cost, 

$ 


$ 0 

$i 

yi 

^2 

I 

6 

8 714 

5. 005X10 5 

0.45 

0. 55 

22 200 

17 920 

150 537 


n 

85 

-900 

5.976x10 s 

.200 

.675 

1 630 

1 029 

7 729 


m 

7 

-13 500 

1. 598x1 0 6 

0.000 

.625 

613 

548 

27 297 


i to m 

98 

-4 452 

5. 392xl0 S 

.550 

.650 

11 520 

6 272 

17 870 



Y-i Yo 

(b) For relation $ = $ Q + $i B max d d equation (17) 


Table 

Number 

of 

coils 

Scaling law 
parameter 

Scaling law 
exponent 

Root 

mean 

square 

error, 

Aq$ 

Mean 

error, 

A x $ 

Mean 

coil 

cost, 

$ 

Ratio of 
scaling law 
exponents, 

6 = y 2 /y x 

$ 0 

$i 

n 


I 

6 

33 840 

2. 426x1 0 4 

1.60 

1.60 

23 040 

18 620 

150 537 

1.00 

n 

85 

1 655 

1.567xlO S 

1.20 

2.00 

2 343 

1 601 

7 729 

1.67 

i to m 

98 

-1 400 

4. 897xl0 4 

1.250 

1.300 

9 581 

5 025 

17 870 

1.04 
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■f _ 'h 

where $ k is the cost of the k coil calculated from the scaling laws, and is the i 

actual cost of the k** 1 coil from the tables. The exponents y^ and y 2 were systemati- 
cally varied to find the values of $q, $j, y^, and y 2 that gave the lowest possible mean 
square error for a particular scaling law. The scaling law of equation (19) gave a far 
worse fit than the other two under all situations, and will not be discussed further in this 
report. 

The values of $ Q , $p y^, and y 2 that give a best fit to equation (18) are shown in 
table IV(a). The four rows of this table correspond to the data of tables I, n, and HI, in- 
dividually and combined. The mean error and the mean square error of the best fit to 
these sets of data are shown, as are the average coil costs for each group. Table IV(b) 
shows the same parameters for equation (17). The computer program would not converge 
to a best -fitting solution for the seven coils of table HI in this case. 

The parameters that give a best fit to equation (18) were used to obtain calculated 
costs for the coils of tables I to HI. These calculated costs have been compared with the 
actual costs of the various classes of superconducting coils in figure 3. The parameters 
from table IV(b) were substituted into equation (17) to obtain calculated costs of the vari- 
ous classes of coils listed in tables I to HI. These calculated costs are compared with 
the actual costs in figure 4. In this figure, the scaling laws would give a perfectly accu- 
rate estimate of the coil cost if the points lay on the diagonal line. In figures 3 and 4, it 
must be remembered that the least -squares curve -fitting procedure used tends to force a 
fit to the larger, more expensive coils rather than the larger number of less expensive 
coils. The dispersion of the calculated costs from the actual costs, even for the best- 
fitting scaling laws, results from the pressure of many marketplace considerations, 
which have little to do with the size of the coil. The magnitude of this dispersion may be 
estimated by comparing the mean error with the mean coil costs in table IV. 

From table IV, it is clear that the mean errors and mean square errors of the scal- 
ing laws of equations (17) and (18) are not significantly different, and either one could be 
employed to provide a rough estimate of the cost of a given coil. In the subsequent dis- 
cussion, equation (17) will be used, since the mean coil diameter is apt to be known at a 
much earlier stage in the design process than the volume of the coil windings, which is 
required in equation (18). 

The parameters of table 1(b) represent the state of the marketplace as of August 1965. 
It is probably not unreasonable to expect that the cost calculated by equation (17) is an 
upper limit to the future cost of superconducting coils, since increased experience and 
improved superconducting materials should tend in the direction of decreased costs. The 
constants $q and would probably be most sensitive to advances in the state of the 
art, but the values of the parameters y^ and y 2 ought to be less sensitive. 
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Calculated cost, dollars 



(c) Niobium zirconium coils in Dewars, Cost of superconducting magnet = 
-1. 35X10 4 + 1. 598X10 6 V°- 625 . 



10 4 10 5 
Actual cost, dollars 

(d) All coils. Cost of superconducting magnet - -4452 + 5.392X10 5 vj?- 65 . 


Figure 3. - Concluded. 


Calculated cost, dollars 



Figure 4 . - Scaling laws calculated from equation (17). Maximum magnetic field on 
axis of mirror machine, B max ; mean diameter of superconducting coil, d. B max 
in Webers per square meter, d in meters. 
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10 3 10 4 10 5 10 6 


Actual cost, dollars 

(c) All coils. Cost of superconducting magnet = -1400 + 4.897X10 4 d 1 * 30 . 

Figure 4. - Concluded. 

Scaling Law for Cost of Conventional, Water-Cooled Copper Coils 

The scaling laws for the cost of conventional water-cooled copper coils have been 
discussed extensively in reference 11. The power required by a conventional coil is given 
in reference 12 as 

pB 2 d 

W = m — watts (22) 

G 2 X(1 + P) 

where p is the conductivity of the conductor, G is the Fabry G factor, \ is the packing 
fraction of the conductor, and P is the ratio of the outer to the inner coil radius. The 
major items of cost in conventional magnet systems are the alternating - to direct -current 
conversion equipment, and the equipment required to supply cooling water. Both are pro- 
portional to the power dissipated by the coils. An approximate scaling law for the cost of 
geometrically similar coils made of the same type of conductor and using the same pack- 
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ing fraction may be written as 


*w = $ ow + $ lw B max d w < 23 ) 

where $ QW is the fixed costs, the constant $ lw is given by 

G^Xd + P) 

and C is the cost of the power supply and water supply per watt. It should be noted that 
equation (24) is of the same form as the scaling law of equation (17) with y^ and 
equal to 2. 0 and 1. 0, respectively. 


GENERAL EXPRESSIONS FOR OPTIMIZED MIRROR MACHINE 

This report is concerned with the optimum design of apparatus for an experiment, 
and not with the plasma properties n Q and V Q , which must be determined by the objec- 
tives (or the results) of the experiment itself. The coil geometry that will be assumed 
for the remainder of this analysis is that shown in figure 5 in which a magnetic bottle is 
formed by two coils of mean diameter d separated by a distance 2Z Q . The magnetic 

field on the axis of the coils is approximated by that of a current loop of radius r =4/2, 

0 


Field-generating 
current loop 


Field-generating 
current loop 


© 


© 



19 



given in reference 13 as 


B z = 


8 (—+ Z 2 


The field B min on the axis of the midplane of the magnetic bottle is 


B . = 
min 


2 ^ 

< 2 o\ 2 / 2 

d 2 + 4Z 2 


and the maximum field B max can be expressed as 


(d 2 *16Z if" 4 


Equations (26) and (27) can then be used to relate the mean coil diameter d to the charac- 
teristic length of the magnetic field Z Q : 


R — 
m 


B min 


d 2 + 4Z 2 


Therefore, 


r =1=Z 0 JL 

2 V R m 


The average magnetic field appearing in equation (13) may be written in terms of the 

mirror ratio R m as 

m mm' max 
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(30) 


B av g + ^max^ ^ B max^ B m + ^ 


H equation (29) is substituted into the figure of merit of equation (12), 

,/l -R_V'*l/„ \2/3 ,_1 

° = n oV* Z ot 


/» - R mV /2 

! ^ 
i 

CO 

V 2 ) 

bn! J 


(31) 


is obtained. When equations (13) and (30) are substituted into equation (14), the adiabatic 
constraint may be written as 


R 


m 


L B max Z o^ R m + 



(32) 


where V Q is taken as the typical particle energy. If equation (29) is substituted into 
equation (17), the scaling law for the cost of a single superconducting coil may then be 
written as 


$ = 


$ 0 + 


SB y l z y 2/2 
*l B max Z o 2 




(33) 


The figure of merit in equation (31) is maximized subject to the constraints of equations 
(32) and (33). 

K B max is eliminated between equations (32) and (33), 


$ = 


$ 0 + 


* /1^2/2-n 

*1 Z Z o 

X'V + D yi 

Am m 



(34) 


K the scaling law exponents are unequal, y^ * can be eliminated between equations 

(31) and (34) and the figure of merit can be written entirely in terms of the mirror ratio 
as 
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Figure of merit, Q 



x T < R m><l ' R n * V] 3ri/(r2 ’ yi> 


(y 2 +2y 1 )/[2(y 2 -yi)] 


After inspection of the terms containing the mirror ratio, it can be verified that, for 
R m = 1.0, the figure of merit S2 is zero. This is expected, since the escape cone occu- 
pies all of velocity space, and no particles are trapped. If the experimental results of 
equation (16) can be extrapolated below R m = 0. 175 to R m « 0, the figure of merit be- 
comes proportional to 


R m-° 


3 rj/ A-(y 2 +2y 1 )/3 


i/(r 2 -ri) 


T <V 



0 Mirror ratio, R m 1 

Figure 6. - Schematic variation of figure of merit with mirror 
ratio. 


The strongest dependence of r on R m is 
such that, unless A is larger than unity, the 
figure of merit will approach zero as the mir- 
ror ratio approaches zero, provided that 
y 2 > y^. From the mean value theorem, it 
is known that SI must have a maximum some- 
where in the range 0 < R m 1. 

When A is less than unity and y 2 > y^, 
the figure of merit has the general form 
shown in curve 1 on figure 6 and possesses a 
very sharp maximum. If y^ > y 2 , however, 
the figure of merit becomes infinite as R m 
approaches zero, as shown on curve 2 of fig- 
ure 6. In the latter case, there exists no op- 
timum mirror ratio, and the figure of merit 
then depends entirely on how low a mirror 
ratio one can afford to pay for. In other 
words, if the coil cost is more strongly de- 
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pendent on magnetic field than on size, there exists no economically optimum mirror ra- 
tio. Inspection of equation (23) shows that this is the case for conventional, water-cooled 
copper coils. Only if the cost depends more strongly on size than on magnetic field will 
an economically optimum mirror ratio exist, and an inspection of table IV (b) shows this 
to be the case for superconducting coils. For reasons discussed on page 30, one would 
expect future progress in the state of the art to result in values of y«/ y l than 

those values shown in table IV(b). 

E equation (35) is differentiated with respect to the mirror ratio and set equal to 
zero, the value of R m for which the figure of merit is a maximum is given by 

1 V 6 - 1 ) R m , 2^ 3 (S + 2) V 6 - 1 ) 3T (37) 

A _ 6(l-R m ) 1 + R m 9 ( 2 ¥3 .rV3) 3t 8R m 

where 6 is defined as 


6 = — (38) 

y l 

Two assumptions are made regarding the variation of confinement time t with mir- 
ror ratio. H it is assumed that the confinement time is independent of mirror ratio, the 
optimum value of the mirror ratio is given by solving equation (37), the last term equal 
to zero: 


*\ C 


^-1 


R. 


6 =^ + 


m 


R 


m 


R 


A 1 + R 


m 


6(1 - R m ) 



>< 


m 


6(1 - R m ) 



(39) 


A theoretical analysis by Lehnert (ref. 14) indicates that this assumption may prove to be 
the case in actual mirror machines, if particles are lost by instabilities or by rapid diffu- 
sion across the field lines in physical space, before they are lost by diffusion into the es- 
cape cone in velocity space. 

If cross -field diffusion in physical space does not occur, the confinement time of an 
ion may be as long as 


T ( R m ) - T co ll P o 


(40) 
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where t cq ^ is a self -collision time that is a function only of the plasma density and tem- w 
perature, and P Q is the loss probability per collision. The confined plasma volume in a 
velocity space of unit radius is just 4 n, and the escape cone occupies a volume 
4ir(l - cos 0 Q ). The fractional volume of velocity space occupied by the escape cone gives 
a rough measure of the loss probability per collision, which is 


P 


o 


1 

1 - cos 6 q 


(41) 


If equations (5) (using B = B min at the midplane) and (41) are substituted into equation 
(40), the confinement time is given by 


T < R m> = 


* coll 

1 - (1 - R ) 
v nr 


1/2 


(42) 


The last term in equation (37) is then 


1 dr 


r dR 


m 


1 

2 


d . R m )l/2 - (1 - R m ) 


(43) 


Under this assumption as to the variation of r with R m , the optimum value of R m is 
given by 



For given values of 6 and A, equation (39) determines an upper bound on the optimum 
mirror ratio, and equation (44) determines a lower bound on the optimum mirror ratio. 
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The true dependence of the confinement time t on R m will undoubtedly lie somewhere 
between that given by equation (42) and the assumption that t is independent of mirror 
ratio. If cross -field ion diffusion is a significant process, the confinement time will be 
somewhat more independent of R m than given by equation (42). 

The optimum values of R ffi given by equations (39) and (44) are of interest to the ap- 
paratus designer, since they are independent of the apparatus dimensions, cost, and 
maximum magnetic field. As pointed out on page 22, however, this optimum, will not 
exist when y^ > that is, when 6 < 1.0. 

The optimum value of R m from either equation (39) or (44) may be used to deter- 
mine the optimum B max and Z Q as a function of the particle energy, particle mass, 
and the available amount of investment capital. Equation (34) can be solved for the opti- 
mum value of Z^, which yields 


Z = 


'» - * o \ 1/h ' 2 '' > ' l) f £ A R m A(R m + « 


0 V *1 


, 6+1 


/ e V/ 2 f 

7 v \^ 3 . r 

6/2 1 

\2mV- / 

\R ) 

( 

V O/ l 

\ m/ 

J 


(45) 


By substituting this into equation (32), the optimum value of B is seen to be 

QlaX 


^mV 
4 ' 0 


A/2 


B max 


6 /( 6 - 1 ) 


$l \i/(r 2 -^i) 


^R^l + R ) 
Am nr 


^ R m/ 


Y /3 - i" 


1/2 


(46) 


Under no circumstances should this optimum value of magnetic field violate the condition 
j3 s 1.0, where 0 is the ratio of plasma to magnetic energy density (to be defined in 
eq. (48)). 

For small coils and some experimental circumstances, the coil cost may not be a 
primary consideration. B: is then of interest to optimize the mirror ratio subject only to 
the adiabatic confinement constraint. K Z Q is eliminated between equations (31) and (32), 
it is found that 


4\& n 0 V 0 « 

7 2ra M 

3 / 2 R 3 t(R )(1 - R ) 1 / 2 
m v m /x m' 

7 2 \ 

2/3 

- 1 

R 3 f 3 
B min € A 

V e ) 

R^/ A (l + R ) 3 
m nr 

\R_ J 

\ m/ 


(47) 
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This equation suggests that the figure of merit may be made as large as desired by de- i 
creasing the magnetic field intensity. However, B max cannot be made arbitrarily small, 
since the magnetic field must be at least large enough to satisfy the condition (ref. 8) for 
the ratio of plasma to magnetic field pressure, which is 


2/i n eV 

a = — — — -< p 

H o K max 

B Z . 
mm 


where 0 is some number less than unity and is determined by stability considera- 

m3x 

tions. If the confinement time is independent of the mirror ratio, it will be given approxi- 
mately by the self -collision time (ref. 15) as 


t ( r ™) * T , 


25. 8 >fn e^m^/ 2 (eV ) 3 / 2 


m' coll 


e n Q In A 


where A is a slowly varying quantity defined by 


A s 1277 / c o V o 


1/2 V e 


If B _ is eliminated between equations (47) and (48) and equation (49) is substituted 
max 

into equation (47), 


103.2 m 2 


3 3/2 5. 5 3/2 , . 

£ A^o e I V lnA 


■fGO 


where the mirror ratio dependence is given by 


[7jl 

RV A ( 1 + R m > 3 LV R » 


It should be noted that, for small values of R m , equation (52) is dominated by the 
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dependence. Therefore, when the coil costs are not a consideration, as small a mirror 
ratio as possible should be used to achieve a high figure ctf merit. 


Possible Extension of General Analysis 

It is tempting to consider the possibility ctf a completely closed-form optimization 
study in which n Q and V Q are also optimized in such a way as to maximize the figure ctf 
merit. In order to include these two variables in the optimization process, however, two 
additional constraint equations are needed, one similar to equation (32) for the adiabatic 
confinement constraint and the other similar to equation (33) for the coil cost constraint. 
Unfortunately, various candidate constraint equations are not as soundly based on obser- 
vational evidence and laboratory experiment as are equations (32) and (33). H they are 
verified by future experimental investigation, various stability criteria might furnish the 
required constraint equations. Equation (48), which states that /3 s 0 , might be used 

as a constraint equation, if and when it becomes known what value of P max is compatible 
with plasma stability and whether and how varies with R m , Z Q , V Q , n Q , etc. 

As soon as valid, experimentally established constraint equations become available, 
it will be possible to extend the present analysis to incorporate an optimization of n Q and 
V Q . The conclusions to be drawn from such an analysis may well differ from the present 
ones, particularly with regard to optimum values off R m and Z Q ; the present state of the 
experimental art, however, does not permit particular values ctf n Q and V Q to be ob- 
tained at will nor does the present understanding ctf plasma theory furnish valid constraint 
equations. As long as this situation lasts, n Q and V c must be left as free parameters 
to be determined by the outcome ctf experimental investigation. 


CHARACTERISTICS OF OPTIMUM MIRROR MACHINE 
UNDER CURRENT CONDITIONS 

It has been shown that if > Vj, there exists an economically optimum mirror ratio 
that is independent of the amount of money expended, the maximum magnetic field, and 
the apparatus dimensions. It is ctf interest to determine the value of this mirror ratio for 
the specific adiabatic and cost parameters given earlier. If the apparatus is designed so 
that the plasma is adiabatic for a Q s tt/2, a generous allowance for the radius of the 
plasma, equations (15a) and (15b) give 


A = 0. 298 


(53a) 


e. = 0.386 


(53b) 





If this value of A is substituted into equation (39), the upper curve of figure 7 of 
6VsR , labeled "r not a function of R , " is obtained. The substitution of this value 
of A into equation (44) yields the curve on figure 7 labeled "t given by equation (42), " 
in which equation (42) gives the dependence of t on R m - 

For a given value of 6, the maximum value of R m is given by the upper curve and 
the minimum by the lower curve. Until experimental evidence on the variation of r as 
a function of R m becomes available, it is not possible to specify the optimum value of 
R m any more closely than to give the range between the two curves on figure 7. 

If the scaling relation of table IV(b) (p. 13) is used for all 98 superconducting coils, 
the scaling law parameters are 


y i =■ 1. 250 


(54a) 
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r 2 = 1.300 


(54b) 


>= 1 ^ 0 = 1.04 

(54c) 

1.25 

$ 0 = -1400 

(54d) 

$ x = 4, 897X10 4 

(54e) 


For the values of 6 shown in table IV(b), the optimum mirror ratio lies in the cross- 
hatched region shown in figure 7. The surprising feature of figure 7 is that the optimum 
mirror ratios are quite high, in the range 


0.95 < R < 1.0 
m 


(55) 


which implies that the economically optimum mirror fields are nearly uniform. 

These high mirror ratios have a large escape cone, and it is unexpected to find that 
such large mirror ratios are in any sense optimum. The basic reason for these hig h 
mirror ratios is that the object of this investigation was not to confine the maximum num- 
ber of particles in a given volume, but rather to confine the maximum number of particles 
of a given energy and density per dollar of investment in the apparatus. The fact that 
these high mirror ratios give a least expensive apparatus should be taken into account in 
the designing of superconducting magnet facilities, especially since current design prac- 
tice for conventional magnet facilities favors mirror ratios of R ^ 0. 5. 

m 

Future progress in reducing the absolute costs of superconducting coils probably will 
not result in the parameter 6 lying outside the range 1.0 ^6 ^2.0. For 6*1.0, the 
cost will be proportional to the magnetic field and diameter raised to the same power, ft 
is difficult to imagine a situation in which the cost of a coil depends so much more strong- 
ly on the dimensions than cm the magnetic field that 6 would be greater than 2. 0, unless 
the coil costs are much smaller than the cost of the cryogenic Dewars which enclose 
them. As can be seen in figure 7, if 1.0 ^ 5 < 2.0, the economically optimum mirror 
ratio will lie in the range 0. 92 ^ R m opt s 1. 00. 

When 6 < 1. 0, no optimum mirror ratio exists, but this situation is not likely to oc- 
cur for superconducting magnets. The maximum magnetic field of a coil is approximately 


B 


~ I - jS 


max 


(56) 


where j is the current density in the superconducting material, and S is the cross- 
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w 

sectional area of the superconducting material. The cost of the coil is approximately 
equal to the cost of the superconducting material, which is 


$ ~ V g ~ dS 


B d 
max 


(57) 


It would be expected, therefore, that 6 = r 2 /r j = 2.0 when j is not a function of the 
magnetic field. In fact, it can be seen in table IV(b) that 6 = 1. 67 for 85 coils of 
niobium - 25-percent zirconium. The value of 5 near unity found for all 98 coils results 
from the influence of the larger, more expensive coils. These coils are designed to oper- 
ate in the region where j is a rapidly decreasing function of B . This dependence 

1 1X1 dA 

of j on B results in values of y. that are nearly equal to y 9 for marginally de- 
signed coils, whose field at the coil windings is so high that j is a rapidly decreasing 
function of the magnetic field. As superconducting materials improve, it should be possi- 
ble to use a material in any given application for which j is not strongly dependent on the 
magnetic field strength. Therefore, values of 6 as large as 2.0 can be expected in the 
future. Even higher values of 6 may occur if the coil Dewars are more expensive than 
the coils. These large values of 6 could come about, because the Dewar costs are a 
function of the coil dimensions and are almost independent of the magnetic field strength. 

It is not difficult to understand now, the appearance of high mirror ratios. The con- 
struction of the figure of merit to represent the total number of confined particles results 

3 

in its proportionality to the volume of the apparatus, that is, to Z . The cost of a coil is 

shown to increase only as a power of Z between 1. 3 and 2. 0, even when B mov is fixed. 

m <ix 

The adiabatic constraint, however, requires that B max z 0 is about constant, so that the 
cost of a barely adiabatic device becomes a weak function of Z Q . The maximization of 
the figure of merit is therefore dominated by the increasing volume obtainable at high 
mirror ratios up to the point at which the reduction in velocity space volume becomes 
counterbalancing. 

It is appropriate to examine the effect on the analysis of the assumptions made con- 
cerning the apparatus geometry. In all previous cases, the cost of a single mirror coil, 
which forms one-half of a magnetic bottle (fig. 2, p. 6), has been discussed. It has also 
been assumed that the magnetic field of the second coil does not contribute to B max on 
the axis of the coil in question. This latter assumption is not a good one at the high mir- 
ror ratios that have been shown to be optimum: it is a conservative assumption, in the 
sense that one will, in all cases, actually get more gauss per dollar than the optimum 
curves imply. Equation (37) shows that the assumption of the magnetic field geometry 
does not significantly affect the optimum mirror ratio. The geometry assumption is em- 
bodied in the third term on the right side of the equation, and this term is small compared 
with the 1/ A term and the term containing t, which dominate equation (37). 
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So far, the implied assumption has 
been made that there will be no impediment 
to using the figure of merit of equation (35). 
It might be the case that the large escape 
cone angles resulting from the large opti- 
mum mirror ratios will result in plasma 
instabilities. These instabilities may make 
it necessary to design apparatus with low 
mirror ratios. It is therefore of interest 
to use equations (51) and (52) to investigate 
how the figure of merit varies when an eco- 
nomic constraint is not applied. Equation 
(51) implies that one can achieve a large 
figure of merit in an adiabatic mirror ma- 
chine by operating at high particle kinetic 
energies and low densities. Equation (52), 
which contains the dependence of the figure 
of merit on mirror ratio, is plotted as a 
function of R m in figure 8, where 
A = 0. 298. There exists no optimum mir- 
ror ratio, since the figure of merit de- 
creases monotonically from infinity at 

II = 0 to zero at Ft =1.0. The depen - 
m m c 

dence of fi on R m is strong, approximately 

1.0 and steeper than exponential for R m <0.5. 

This analysis shows the sharp contrast between the case in which economic factors 

are an important design consideration and the case in which they are not. When y^ > y 2 

and/ or economic constraints are not important, figure 8 makes it clear that one should go 

to the smallest possible mirror ratio. On the other hand, if y 2 > Tp as is the case 

with superconducting coils, and if cost is significant, the maximum figure of merit occurs 

at a high mirror ratio 0. 95 ^ R , ^ 1. 00. 

m, opt 

Since there exist some circumstances in which one will not operate at the optimum 
mirror ratio, it is of some interest to use equation (34) to estimate the cost of a coil of 
given mirror ratio used to confine particles of energy V . If the parameters of equations 
(53) and (54) are substituted in equation (34), 


exponential in the range 0 . 5 R 


m 


-1400 + 14. 9(f V o )°- 625 


,0.05 


Rt 2 (UR m ) U5 ° 

m m' 



-0.65 


(58) 
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where £ is the ion mass in atomic mass units. The dimension Z Q is of the order unity v 
in the mks system of units adopted, and the dependence is so weak that one can define 
~ j o, Equation (58) is plotted on figure 9, which shows the cost of an axisym- 
metric superconducting solenoid required to confine adiabatically particles of mass £ 

(amu) and of energy V o (eV) in a mirror field of the indicated mirror ratios. It is obvi- 
ous at once that it is very expensive to confine energetic particles adiabatically in mir- 
rors with low mirror ratios. For example, in order to confine 50 keV deuterons 
(V Q = 5X10^, £ = 2. 0) in a 2 to 1 mirror ratio field (R m = 0. 5), approximately $170 000 
must be spent on a single mirror coil. 

The information presented in figure 9 can be presented a little differently by re- 
writing equation (56) as 



Figure 9. - Cost of axisymmetric superconducting solenoid Figure 10. - Critical nonadiabatic energy as function of mir- 

as function of mirror ratio. ror ratio and cost 
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( 59 ) 



The particle mass -energy product £V Q is platted as a function of mirror ratio for sev- 
eral values of capital investment cm figure 10. to order that a particle of a given £V 
confined by a given R m be adiabatically confined, the amount invested must be at least 
as much as indicated cm the curves. Again, it is clear that adiabatic confinement at low 
mirror ratios is very expensive, to figure 9, adiabatic trapping will occur if the amount 
of capital invested is equal to or greater than that shown on the energy curves; in fig- 
ure 10, adiabatic trapping will occur only if the particle energy is less than that indicated 
cm the capital investment curves. 


CONCLUSIONS 

A figure of merit for controlled fusion research apparatus has been proposed, the 
maximization of which is consistent with the needs, objectives, and capabilities of current 
controlled fusion research. The proposed figure of merit is proportional to the apparatus 
volume, the degree of isotropy in velocity space, and the charged particle density, en- 
ergy, and confinement time. 

A scaling law for the cost of superconducting solenoids has been derived by a least- 
squares curve fitting to the actual cost, dimensions, and maximum magnetic field of 
98 solenoids that have been sold commercially. This scaling law may be used for pre- 
liminary design studies and refined on the receipt of actual bids for a specific facility. 

The figure of merit has been maximized subject to an economic constraint, which 
was provided by the scaling law for coil cost, and an adiabatic confinement constraint, 
which was provided by the results of a recent experimental investigation. E was found 
that, when superconducting coils are used, there exists an economically optimum mirror 
ratio that is independent of the dimensions and magnetic field of the mirror coil and of 
the amount of capital investment in the superconducting coils. This optimum mirror ra- 
tio permits the adiabatic confinement of the maximum total number of ions of a given en- 
ergy and equivalent isotropic density per dollar of capital investment in the magnet facil- 
ity. to the optimization process, it was assumed that the ion energy and equivalent iso- 
tropic density were not functions of the mirror ratio, of B„ or of Z . It was also 
assumed (in keeping with the current state of the art) that the ion energy and density were 
dependent variables whose numerical values depend on the outcome of an experiment, 
and not independent variables, whose values can be adjusted at will. 
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The economically optimum mirror ratio was calculated for adiabatic confinement in * 

an axisymmetric mirror field generated by superconducting coils. It was found that the 

economically optimum mirror ratio was quite high, in the range 0. 95 ^ B . / B < 

mirr max 

1. 00. This optimum is a sharp maximum in the figure of merit. Thus, if stability con- 
siderations and other factors permit, designers of superconducting magnet facilities for 
controlled fusion research should aim for much higher mirror ratios than have been cus- 
tomary in existing conventional magnet facilities. 

The large anisotropy in velocity space associated with the near -unity optimum mirror 
ratios may cause such severe instabilities that lower mirror ratios (and, hence, greater 
isotropy in velocity space) will be mandatory. It may also be desirable to use off -optimum 

values of R , B m , v) or for other reasons. The two constraint equations are writ- 
nr max’ o ^ 

ten in such a way as to give the cost of a superconducting mirror coil as a function of the 
confined particle mass and energy and of the mirror ratio (not necessarily optimum) of 
the apparatus. R is shown that the cost of going to adiabatic mirror ratios less than 0. 50 
is quite high for ions of thermonuclear temperatures and will remain so even if the cost 
of superconducting coils is substantially reduced. 

When cost is not a major consideration in designing a magnet facility, the facility 
should be designed with as low a mirror ratio as is consistent with adiabatic confinement. 
When cost is a major consideration, superconducting coils should be designed with a mir- 
ror ratio of about B min /B max = 0. 95. This conclusion is virtually independent of the 
cost of superconducting coils at any given time. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 31, 1966. 



APPENDIX - SYMBOLS 


[Unless otherwise noted, the rationalized mks system of units has been used 
throughout this paper.] 


A 

parameter given by equation (15) 

P 

ratio of outer to inner coil radius 

B 

2 

magnetic field strength, W/m 

P o 

defined by equation (41) 

B av 

average magnetic field on axis, 
»/2 < B mm + B max>' w /“ 2 

R m 

mirror ratio B . /B _ on 

mur max 

apparatus axis 

B max 

maximum magnetic field on axis 

A 

r 

radial distance in configuration 


of mirror machine, W/m * 


space 

B - 
mm 

minimum magnetic field on axis 

r c 

mean radius of mirror coils 


of mirror machine, W/m^ 

r p 

mean radius of plasma 

C 

cost of power supply, dollars 

t 

time 


per watt 


d 

mean diameter of superconduct- 

V 

particle energy, eV 


ing coil, m 

V o 

particle energy on axis where 


mean diameter of conventional 


B = B . , eV 
min’ 


coil, m 

V s 

volume of superconducting wind- 

e 

electric charge, C 


ings 

f 

velocity distribution function of 

V 

particle velocity 


ions 

W 

power required for conventional 

G 

Fabry G factor (see ref. 12) 


coils 

I 

current, A 

Z 

axial distance in configuration 

3 

2 

current density, A/m 


space 

K 

ratio of mean coil to mean 

Z o 

axial distance between B - 

min 

and B T¥ , ov 
max 


plasma radius, r c /r p 


a 

particle mass, amu 

a 

dimensionless radial coordinate 
defined by equation (9b) 

m 

particle mass 

% 

dimensionless radial coordinate 

n 

ion density (see p. 5) 

defined by equation (9c) 

n o 

ion density on axis where 

0 

ratio of plasma to magnetic pres- 


B = B min (see SUMMARY) 


sure, equation (48) 
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^max 

defined by equation (48) 

X 

packing factor of conductor 

r l 

scaling law exponent in equations 
(17) to (19) 

** o 

permeability of free space, 
4 ttX10 -7 H/M 


scaling law exponent in equations 
(17) to (19) 


geometry factor defined by equa- 
tion (11) 

Aj$ 

mean error of scaling law, de- 

P 

conductivity of conductor 


fined by equation (20) 

T 

characteristic confinement time 

A 0 $ 

mean square error of scaling 


of ions 


law, defined by equation (21) 

<P 

azimuthal angle of spherical co- 

6 

parameter defined by equa- 
tion (38) 


ordinate system in velocity 
space 

e 

adiabatic parameter defined by 

P 

volume of plasma 


equation (13) 

ft 

figure of merit (defined by eq. 

e A 

parameter defined by equation 


(D) 


(15b) 

$ 

cost of superconducting magnet 

e o 

permittivity of free space, 
8. 854X10" 12 F/m 

*A 

actual cost of superconducting 
magnet 

V 

dimensionless axial coordinate 
defined in equation (9a) 

$o 

scaling law parameter in equa- 
tions (17) to (19) 

e 

polar angle of spherical coordi- 
nate system in velocity space 

K 

$i 

cost of conventional magnet 
scaling law parameter in equa- 

9 o 

o' 

angle of adiabatic escape cone in 
velocity space given by equa- 
tion (5) 

azimuthal angle in configuration 
space 


tions (17) to (19) 


fe 
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